RADIO SCIENCE Journal of Research NBS/USNC-URSI 
Vol. 69D, No. 3, March 1965 

Propagation in Nonuniform Gyrotropic Media ^ 

Stanley H. Gross 
Airborne Instruments Laboratory Division of Cutler-Hammer, Inc., Long Island, N.Y. 

Leopold B. Felsen 

Department of Electrophysics, Polytechnic Institute of Brooklyn, Brooklyn, N.Y. 

(Received September 8, 1964; revised September 24, November 4, 1964) 

With the exception of the stratified layer approximation, the problem of wave propagation due to 
localized sources in nonuniform magneto-ionic media has received little attention. In previous work, 
the coupling between natural modes of the homogeneous gyrotropic medium has restricted the useful 
range of the analysis to the slowly varying case. In the present investigation, an E and H mode repre- 
sentation is used which results in a type of coupling independent of the rate of variation in the medium, 
thereby permitting the study of regions with strong inhomogeneities along the gyrotropic axis. 

The coupling between the E and H modes now arises due to the deviation of the gyromagnetic 
parameter Y=(Oclco from zero or infinite values. The uncoupled problems are therefore associated 
with the nongyrotropic {Y=0 or Y—^) but nonuniform cases, and the procedure converges rapidly 
when Y is large or small, even for rapid spatial variations of the medium parameters. 

An iterative expansion is employed to deduce the gyrotropic corrections, and the procedure is 
phrased systematically in transmission line form. The analysis is applied to the evaluation of the far 
field radiated by an electric current element exterior to a semi-infinite, nonuniform plasma region. 

1. Introduction 

A fundamental difficulty in the analysis of wave propagation in a continuously stratified mag- 
neto-ionic medium is due to the coupling among the four wave types (upgoing and downgoing, 
ordinary and extraordinary) which may propagate independently when the medium is homogeneous. 
Four coupled first-order differential equations may be derived for the specification of the wave 
functions [Budden, 1961; Clemmow and Heading, 1954] and while a formal exact solution may be 
found under special conditions [H. Keller, 1953], the determination of explicit results has generally 
been based on perturbation techniques which converge rapidly when the coupling is small. Such 
approximation procedures have utilized as unperturbed solutions the magneto-ionic waves in a 
homogeneous region, with the coupling produced by the variation of the medium parameters 
along the axis of stratification. While the coupling is minimized when the variation is small, 
thereby permitting an approximate analysis of the four equations for the wave functions [Clemmow 
and Heading, 1954], there exist regions even under these conditions wherein two or more wave types 
may interact strongly (for example, near reflection points). A modification may then be employed 
by which the problem is phrased generally in terms of two coupled second-order equations 
[Hougardy, 1961]; the two solutions of each equation (for example, upgoing and downgoing waves 
of one type) may then t)e strongly coupled, whereas, interaction between the two pairs is weak. 
The analysis in all of these considerations is built around the unperturbed solutions for the homo- 
geneous case, thereby making slow variation an essential requirement. 



'The material for this paper has been taken from the dissertation of one of the authors, S. H. Gross, submitted to the Faculty of the Polytechnic Institute of 
Brooklyn in partial fulfillment of the requirements for the degree of Doctor of Philosophy, 1964. 
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An alternative viewpoint is employed in the present treatment where the unperturbed situation 
is one for which the magneto-ionic dielectric tensor has a diagonal representation. In a coordinate 
system oriented so that the z-axis coincides with the direction of the d-c magnetic field B, this 
occurs when B = ^ (uniaxial case) and, of course, also when B = (isotropic case). If the uniaxial 
medium parameters vary only (but arbitrarily) along 2, four uncoupled wave solutions may be 
derived; the coupling in this instance is due to finite, rather than infinite, values of B (or of the 
gyrofrequency coc). The present formulation is therefore useful for the study of the influence 
of strong gradients in the constitutive parameters provided, however, that coc is large in comparison 
to the applied frequency co and the plasma frequency cop. (The constraint on cop relative to a>c is 
actually less severe, see sec. 2.3b.) While this regime may exhibit strong anisotropic effects, 
analogous considerations apply as well to the almost isotropic range where coc is small. 

The analysis proceeds from a guided wave viewpoint by which the field is decomposed into 
E and H modes with respect to the direction of B. The modal amplitudes are found to satisfy 
conventional transmission line equations which are coupled by the off-diagonal term in the dielectric 
tensor. A solution of these equations is achieved by an iterative scheme which converges rapidly 
when coc is very large or small. The procedure is illustrated for the case of a transverse electric 
current dipole embedded in a magneto-ionic medium with z-dependent electron density, with field 
expressions derived for the unperturbed situation (coc = ^) as well as the first gyrotropic correction, 
O(l/coc)- Explicit results are summarized for the case of a half-space with constant density gradient, 
with the source located in the exterior region. 

A variety of other applications for Y > \ and Y < \ have also been studied in detail and will 
be submitted for publication in separate form. It may be mentioned that although the case 
Y <\ pertains more directly to the ionosphere than Y > 1, the latter has been chosen for illustration, 
since even the unperturbed solution (F=oo) may exhibit strong anisotropic effects. 



2. Formal Solution 

As mentioned earlier, the direction of nonuniformity and the steady magnetic field are taken 
along the z-axis of a right-handed Cartesian coordinate system, so that the configuration in the 
x-y plane is uniform and infinite in extent. The relative permittivity tensor "e relating the electric 
displacement vector and the electric field intensity vector in the nonuniform, gyrotropic medium 
may be expressed as foUows [Budden, 1961]: 
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(1) 

where the double bar above the symbol denotes a tensor, the expression on the right is its matrix 
representation and the elements of the matrix ei, €2, and €3 are permitted to be functions of z. 
The permeability of the medium is taken as that of free space. Throughout this paper the implied 
time dependence is taken as e^^^^^ where co is the radian frequency of the source, and the steady 
magnetic field is taken along the negative z-axis. 

The quantities ei, €2, and €3 for a lossless gyrotropic plasma medium have the form [Budden, 
1961] 
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^here: 

X = a)p/(jt)^ 

Y=(x)clo) 



cop = plasma frequency = \/ 

' ^' eom 



o)c = gyromagnetic frequency = 



^ = steady magnetic field, 



eB 



€o is the permittivity of free space, e and m are the electronic charge and mass and A^ is the electron 
density of the medium. 



2.1. Reduction of the Maxwell Field Equations 

It has been shown elsewhere [Arbel and Felsen, 1963] that the Maxwell equations for the 
transverse field components may be written as follows: 



-^ = -iojeoil^]T,VtVt)'ZoXEt^zoXj, 
dz k^ 



- ^ = iwii^Cl ^ + 77- V, V,) • zo xHt- zo X M, 
dz Ic^e^i 



(3) 



n 

where the subscript t is used to designate vector components transverse to z, V/ = V — Zq— , 

k^ = (o^lJio€o, E and H are the electric and magnetic field intensity vectors, /xo is the permeability 
of free space, 1^ is the transverse unit dyadic, Zo is the unit vector along z, and it is the transverse 
tensor givey by: 



_. e\ — ie 



;) 



Furthermore, 



- - 1 

Jt = Jt-j—-VtM,Xzo, 

lOJfJio 



- - 1 

Mt = M(^- VtJzXzo, 

l(0€o€',i 



are equivalent transverse source distributions whose values are expressed in terms of those of 
the specified electric and magnetic vector current densities J and M, respectively. The z com- 
ponents of the field vectors, Ez and Hz, may be expressed in terms of the transverse field com- 
ponents and the z components of the applied current densities [Arbel and Felsen, 1963]. 

To solve_[3) in a_transversely unbounded region, one seeks a representation for the transverse 
field vectors Et and Ht in terms of a complete set of vector eigenfunctions. If these eigenfunctions 
belong to a ''characteristic" set for the configuration in question, they are solutions of the source- 
free field equations, satisfy simple orthogonahty conditions, and their z-dependence in a homo- 
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geneous medium is of the form exp( ± iKnz), where Kn is the longitudinal wave number and the 
subscript distinguishes the nih mode. Such a set of functions has been found by Arbel and 
Felsen [1963] and utilized by them for the representation of fields in piecewise homogeneous 
regions; the resulting "ordinary" and "extraordinary" modes, familiar from conventional magneto- 
ionic theory, are uncoupled in each homogeneous section but couple at the boundaries. If the 
eigenfunctions turn out to be independent of the properties of the medium in the transverse plane, 
then they lead to an uncoupled description even in the presence of continuous stratification. This 
is not the case for the characteristic modes and leads to the coupling difficulties mentioned in the 
introduction. 

By the alternative approach explored in this paper, the mode set is chosen so as to satisfy the 
source-free field equations when the dielectric tensor is diagonal (uniaxial or isotropic case). In 
this instance, the resulting eigenfunctions are indeed independent of the transverse medium para- 
meters (they are in fact identical to the transverse E and H mode functions for isotropic regions) 
and thus lead to an uncoupled representation for arbitrary ei(z) and 63(2) when €2 — 0. 

We therefore assume the field representation [Felsen, 1964, sec. 2.2] 

n n I 

> (4) 

Ht = xnizihup) + 2/;'U)A/;(p), J 

n n 

where em and hm are the transverse E and H mode functions [see C^) and (15)], the F's and /'s are 
voltage and current coefficients (yet to be determined), htn = ZoXetn, p is the transverse position 
vector, and the single and double primes refer to E modes and H modes, respectively. The 
summation is over all values of the modal index n (a double index), which is wholly or partially 
continuous for the present case of an unbounded cross section. For example, if the cross section 
is described in cylindrical coordinates the eigenfunctions and eigenvalues take the form given 
in (15) and the first paragraph of section 2.3a. 

On inserting (4) in place of the transverse field components in (3), and using the known orthog- 
onality properties of the eigenvectors and the equivalence e/^ = h"^, ^tn~~^tn ^^^ modes in an un- 
bounded transverse cross section [Felsen, 1964], one obtains the following two sets of differential 
equations for the modal amplitudes: 
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where Kn and Zen are the modal propagation constant and characteristic impedance, respectively, 
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and ktn is the transverse eigenvalue corresponding to the modal index n. (In a cyhndrical coordi- 
nate representation, ktn is equal to the radial eigenvalue kt.) The quantities i^, v'^, i"^, and v" are 
the modal voltage and current sources which are obtained from the current densities J and M as 
follows: 



i;= // dSJrel^X' 



fS dSJre';:, 
s 



v',= ff dSMfh'*,v:= IS dSMrK*, 



f^tn f 



(7) 



where the integrals are taken over that portion of the transverse cross section which is occupied 
by the sources. (/)« denotes a normalized scalar eigenfunction (see beginning of sec. 2.3a). 

Equation (5) contains the differential equations for the E mode voltage and current coeffi- 
cients, whereas (6) is for the H mode quantities. Each set contains a term which cross couples 
the modes, the coupling coefficient being the same in both cases, namely 0)6062. Consequently, 
the cross coupling is proportional to the degree of gyrotropy of the medium as represented by 62 
and exists whenever this parameter is nonvanishing. When 62 = 0, (5) and (6) reduce to the uni- 
axial forms [Felsen, 1964] and are uncoupled. For the isotropic medium, 61 = 63 and 62 = 0, in 
which instance the modes are uncoupled and k^^ = k^, as expected. Since 61, 62, and 63 are func- 
tions of 2, (5) and (6) are differential equations with variable coefficients. They may be inter- 
preted in network form, with V,/ and In representing the voltage and current on an E mode non- 
uniform transmission line; V,/' and 1/ denote the corresponding H mode quantities. In such a 
representation the lines are coupled everywhere along their lengths, and are excited by the genera- 
tors, in\ Vn , in\ and Vn" . 

Since 62 in (2) contains X and Y but not their derivatives, cross coupling in the E and H mode 
representation does not depend explicitly on the rate of variation of the medium. This result is 
unlike that obtained from the formulation in terms of ordinary and extraordinary modes in which 
the rate of variation enters directly into the coupHng terms [Fosterling, 1942; Clemmow and 
Heading, 1954]. Consequently, the E and H mode coupling coefficient is not necessarily large 
for large gradients, and the utility of the solution need not be restricted to the slowly varying case. 

Upon solving (5) and (6) simultaneously for the modal voltages and currents (subject to as yet 
unspecified boundary conditions) and then placing these quantities into (4) we obtain in principle, 
the formal solution for the transverse fields. The z components are then derived via [Arbel and 
Felsen, 1963] : 



E, = - 



H,= 



1 



tct)6o63 n 

1 



^ILVt-e'tn^^ 



I0)eo€3 



icojjio n 



^VlVt'^tn^^ 



M, 



KOjJLo 



(8) 



Because of the complication introduced by the cross coupling terms in (5) and (6), it is in general 
quite difficult to provide explicit expressions for the field components in a nonuniform, gyrotropic 
medium. However, by making 62 small enough, the equations are amenable to treatment by an 
iterative process, as will be demonstrated. It may also be noted [Gross, 1964] that it is possible 
to transform (5) and (6) so that the resultant modes are uncoupled and independent when the 
medium is assumed uniform but gyrotropic. The solution in this case is the same as that found by 
Arbel and Felsen [1963]. 
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2.2. Iterative Expansion 

From (2), it is seen that €2 will be small if F <^ 1 or Y> 1. The expansion will be developed for 
the latter case only, but may be carried out in an analogous fashion for the former case [Gross, 
1964]. 

Assuming that F is a very large but finite constant (i.e., the inhomogeneity rests solely in the 
parameter X), we represent the modal voltage and current coefficients in series of powers of l/F 
as follows (the modal indices n will be omitted for convenience): 

^^ -^^y^' ^= ^ W (9) 

m=0 -«■ m = JL 

and a similar expansion is performed for ei and €2 (63 is independent of Y). On inserting these 
expansions into (5) and (6) and equating terms with equal powers of 1/F, we obtain a set of dif- 
ferential equations for each iterative order. The zeroth order equations are those for the uniaxial 
voltage and current coefficients, 

dV" dl" i ' 

whereas, the mth order gyrotropic corrections to the uniaxial coefficients are specified by (note: 
quantities with negative index vanish identically), 

dV;,, . k^{l-X)-kf 
dz ' a>eo(l-Z) '" ' 

df 



—7^ — iw/Ao/m = 

dn, 



(11) 



dz (jofjio 

It is seen that the true sources appear only in the zeroth order solution and that the source for the 
mih order coefficients is provided by the voltages of lower order. Since the homogeneous forms 
of (10) and (11) are identical, the solutions for any order may be derived from the same Green's 
function, and the E and H mode equations are no longer coupled. The complete solution may 
therefore be obtained, in principle, by an iterative procedure once the uniaxial result has been 
determined. It may be mentioned that the H mode coefficients are independent of the medium, 
in consequence of the perturbation about the uniaxial case which does not influence the H mode 
fields [Felsen, 1964]. 



2.3 Excitation by a Transverse Electric Current Element 

a. Voltage and Current Coefficients 

To illustrate the form of the solution obtained from the preceding analysis, we consider the 
case of a transverse electric current element located at z = z' on the 2-axis and oriented parallel 
to the X-axis. Cylindrical coordinates (p, i//, z) will be used, where p is the radial distance from 
the z-axis in the transverse plane, and i// is the angle between the radius vector p and the x-Sixis. 
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In this coordinate system the modal index n is described by the two quantities kt and v which 
represents the eigenvalues in the angular and radial domains, respectively, and the normalized 
scalar mode function is given by 

(t)n = {ktl27TVl^ exp {iv\lf)J^{ktp), ^' = 0, ±1, ±2 . . ., < A;, < oo 

(Felsen, 1964, part 2). For a tranverse current source, located as mentioned above, the only non- 
vanishing modes are those for which v = ±.\. Furthermore, a transverse electric current element 
gives rise to vii =Vn =0 for all modes, as may be seen from (7). Then with J = xo8{z — z')8{p), 
where Xo is a unit vector in the x-direction and 8{a) is the delta function, we find from (7) (see 
also (15)): 



i;;-a;;8(z-2')=T|^8(2-z') 



i';^^a:fi(z-z')^^^8{z-z') 



(12) 



where it is understood that n now signifies the double subscript (kt, ±1). al^ and aj| denote the 
strengths of the current generators located on the transmission lines at z' . 

For the zero order E mode coefficients [F^^ and 7^^^, with m = in (9)], 

f;„(z,2')=-<z;.,(z,z') 

(13) 
i:,„{z,z')=-a:n',(z,z') 

where Z^^(z, z') and r^J(2, z') are voltage and current Green's functions [Felsen, 1964, sec. 2.3] 
for the uncoupled E mode uniaxial transmission lines excited by a current generator of unit strength 
at z' . The same expressions apply for H modes, with double primes replacing the single primes. 
Since the Green's functions depend on kt but not on v [see (22) and (23), with /c, F, and Z given 
in (6)], the latter index has been omitted in the suffix. In the first order equations, the zero order 
voltage coefficients V^^ and FJJ^ act as continuously distributed current generators, as may be 
seen from (11). Thus, the total voltage and current coefficients for the modes n = {kt,±il) are to 

o(i/y), 

I I k, a.€o J "^ j(^")^^^(^^ z")Zl{z!\ z')dz" + 0( 1/P) 
n[{z,z')-^yl^^jy{z'')Tl'lz,z'')Ziy,z')dz'' + 0{llY^) 



2 V27r 




(14) 



^ Zllz, z') +i yj^ ^ jy{z")Zl{z, z")Z;/z", z')dz" + 0(1/P) 



Where multiple signs occur, the upper sign goes with v = -\-l, and the lower with v = —l. One 
observes that the transverse electric source excites gyrotropic corrections of the order 1/F in both 
the E and H mode coefficients. This behavior occurs also with a transverse magnetic current 
source, whereas for longitudinal sources the gyrotropic effects are of higher order [Gross, 1964], 
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For a longitudinal electric current element the first E mode gyrotropic correction is in the l/F^ 
term, while the total H mode coefficients are of 0(1/7). 

For the nearly isotropic medium (Y « 1), similar results can be derived readily. 

b. Field Components 
In a cylindrical geometry the vector eigenfunctions are given by 

4,^, (p, i//)=X';,,, (p, i//)--^ e'-^{po |(A-i(M-A+i(M)+^o|^y#^^^ 

where Jp is the Bessel function, and po and i//o are unit vectors in the p and i// directions, respec- 
tively. The field components are obtained from (4), (8), and (14), with y]"^" ^ I dktkt; only 

the v = ±^ terms contribute when the z-axis passes through the source. The following expressions 
result for the field components in a cylindrical coordinate representation: 



^P* = ^ /o *^' { si°n J [+-/«(M(^;,U, z') + Zl{z, z')) + Uhp){Z;^{z, z')-Zl^{z, z' 



))] 



0)60 sin \\} 
Y cos i// 



J,{k,p) J" dz"X{z") {Z',,(z, z")Zl,{z", z') + Zl,{z, 2")Z;/2", z')) 
-jAktp) \" dz"X{z"){Z'^^{z, z")Zl^{z!\ z')-Zllz, z")Z',lz", z') 

J -co 



+ 0(1/P) (16) 



where the upper algebraic signs and trigonometric functions refer to Ep, and the lower symbols 
are appropriate to £'^. 
Also, 

E^^ ^TTcJeoe, j^dhkUiikp)^mz, z') cos ilj^^'jyz"X{z")rij^z, z")Zl,{z\ z')}+0(l/F2). (17) 
Similarly, 

^p^=-^ If *^^ { JI3" J[jo(M(n;(^, z')+7i;(^, z'))±Mktp){Ti]{z, z')-v^^{z, z'))] 



_ ajeo cos i// 
F sin i// 



yo(i,p) ^^dz"X(z"){TUz, z'')Zllz", z') + T;Iz,z")Z '„(z", z')) 
±j2{hp) jyz"X{z") {n,{z, 2!')Zl,{z\ z') - T';,(z, z")Z;,,(z", z')] } + 0(1/P). (18) 

Finally, 

//. = ^^ ^ykckUAktp) {zr,(z, z')-^ ^'^dz"X{z")Zl(z, z")Z'4z", z')} +0(1/F). (19) 

The zero order expressions (in 1/F) in these equations yield the uniaxial field components in 
terms of uncoupled E and H modes, whereas the first order terms give the first gyrotropic correc- 
tions. The solutions for the field components in (16) to (19) must still be regarded as formal since 
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it is necessary to sul)stitute for the Green's functions and carry out the integration. This will be 
demonstrated below in a simple example. More complex problems have also been examined 
[Gross, 1964] and will be submitted separately. 

Since the iterative solution involves a power series in l/F, there may be some question con- 
cerning its convergence. In view of the complicated nature of the coefficients in the series, a 
general statement is difficult and we must trust the expectation that reliable results may be ob- 
tained with a few terms when Y is sufficiently large. The following observation is helpful in this 
connection: since each integral in the iterated solution contains the factor A^(z) (the mth order term 
involves an m-iold integration), the mlh order term may behave according to (Xm/V)"^ where Xm 
denotes the maximum value of ^(z) in the integration interval. Thus, it is to be anticipated that 
the expansion is useful when (XmlY) «l. The restriction (XmJY) < I is well motivated since 
this range, but not the complementary one, may be regarded as a perturbation about the uniaxial 
case. In certain configurations, such as the semi-infinite medium with monotonically increasing 
electron density, X may reach values in excess of Y for distant regions along the z-axis. If slight 
losses are included in these cases, it can be shown [Gross, 1964] that these regions will not limit 
the convergence. Instead the rate of variation of the medium parameters over a free space 
wavelength must be made small compared to Y; this condition appears in (44) as ^ X << 27tY. 



2.4. Construction of the Modal Green's Functions 

As noted in sections 2.2 and 2.3, the iterative solution may be built up entirely from the Green's 
functions for the transmission line equations (10) or (11), the functions being identical since these 
equations differ only in their source terms. The various Green's functions have a direct physical 
interpretation: Yn(z, z') and Tn\z^ z') denote the current and voltage, respectively, observed at 
z due to excitation by a unit strength voltage point source at z', while Zniz^z') and Tn{z^z') denote 
the voltage and current, respectively, observed at z due to excitation by a unit strength current 
source at z' . While the transverse electric dipole element discussed above does not give rise to 
a voltage source, thereby requiring the use of Zn(z, z') and r,/(z, z') only, all of the Green's func- 
tions are needed for a description of the field with more general excitation. Since the voltage and 
current at a source-free point are connected by the transmission fine equations, 

= IKnZcnIniz, Z ), ^ = iKn/r/t^wl-Z, Z ) (20) 

and in view of the reciprocity relations [Felsen and Marcuvitz, 1956, 1959] 

T^iz, z') = - JZ{z\ z), Yniz, z') = Y„{z\ z), Z„(z, z') = Z,(z' , z), (21) 

all of the Green's functions may be derived from the knowledge of any particular one. In (20), 
Kn, Zen, and Ycn=VZcn represent the previously defined propagation constant, characteristic im- 
pedance and characteristic admittance, respectively, all of which may be z-dependent (the iden- 
tification of these quantities in (10) or (11) is evident). It is usually convenient to determine 
either Yn{z, z') or Zn{z^ z'), with the former and latter preferable for the E and H mode problems, 
respectively. In terms of Yn{z, z'), 

^"^^' ^'^ K„{z)Yrn{z)K,.{z')Yeniz') dzdz' ^"^^' ^'^ 

with dual relations applicable when Z„(2, 2') is given at the outset. 
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Figure 1. Nonuniform section of transmission line. 



The general network problem for the determination of Y{z, z') is shown in figure 1 (the modal 
subscript n is omitted henceforth); in mathematical terms, Y{z, z') satisfies the second-order dif- 
ferential equation obtained after elimination of P from the transmission line equations: 



d 1 d 



dz k{z)Yc{z) dz 



- k{z)Zc{z)\y{z, z') =ib{z- z') , (23) 



The nonuniform transmission line described by the parameters k{z)^ Zc{z) is assumed to be termi- 
nated at zl and zr by impedances Ztl and Ztr^ respectively, which prescription specifies the bound- 
ary conditions at the z-termini of the region covered by (23).^ The construction of the Green's 
function is conventional: we assume that the two linearly independent solutions of the homogeneous 
equation (23) are superposed in such a manner as to yield two (current) functions Il{z) and Ir(z) 
which satisfy the boundary conditions at zl and zr, respectively. Then 

y, ,._ Il{z<)Ir{z>) . . diR J dh 

where 2< equals z when z <z' and z' when z > z\ with the converse applying to z>. The denomi- 
nator in (24) (the Wronskian) is a constant. Upon introducing the impedances, 

77 \^ Vl{z) ^ dJL{z)ldz -^, . ^ dlR{z)ldz .^r^ 

^^^^ h{z) iK{z)Yc(z)h{z)' ^^^^ iK{z)Yc{z)lR{z)' ^^^ 



Y^^ j^^ [hiz<)llL{z')] [Ir{z^)IIr{z')] 



seen to the left and right, respectively, of the point z, one may write (24) in the normalized form, 

')] [Ir 

Ziz') 

(26) 
where 

Z(z')=Z(z')+^(z'). 

The detailed structure of II and Ir may be exhibited either in a standing wave FFelsen and 
Marcuvitz, 1959] or a traveling wave representation. For the nonuniform transmission line prob- 
lem in question, the latter is frequently more convenient and provides an insight into the physical 
mechanism of reflection and refraction in a progressing wave. We shall therefore break the total 
currents into two parts 

h=ii+ii,h=n+i], (27) 



^ Evidently, no mode coupling is assumed to take place at the boundary. This obtains when the regions z < Zi and z> zr are occupied by an isotropic or by 
different anisotropic medium (with optic axis along z), as long as the analysis is performed by the £ 



and H mode decomposition. 
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with superscripts -f and — denoting waves traveling in the (-^z) and (—z) directions, respectively. 
While these distinctions become obscured in an inhomogeneous region where continuous reflec- 
tions take place, they may be made precise by the requirement that the + wave satisfy the boundary 
conditions at z— >«» (radiation condition) if the medium is infinite in extent, with an analogous 
requirement imposed on the — wave at 2 ^ — oo. Thus, 

^. . dl-{z)/dz ^. dl^{z)ldz 

^'^'^=~iKYJ-{z) ' ^'^'^^JkYJ^) (28) 

represent the input impedances to the unbounded medium (matched case) seen to the left and right 
of 2, respectively, and 

*- . _ Itiz) _ t{z)-Z'{z) ^, , _ /«(2) _ Uz) -tjz) 

i. (z) ———— — — , lu;— —— ~ — 

IZ(z) Zi{z)+Z{z) I^iz) Zi(z)+Z{z) (29) 

are the (current) reflection coefficients. Tiie reflection coefficients at any two points may be re- 
lated by the identity 



^ ^ Ijiz') Il(z) 

r{z)= nz')-fjj^jzj^, 

with a similar expression for r{z), and it is then not difficult to show that 

j,(^ ^,) _ [/,:(z<) + n2,.)/M2<)] [h(z>) + r(zSi',(z>)] 

Zi(z')[/,^(z')/,:(2') -?(2«)r"(2,.)/,7(2')/,t(z')] 



(30) 



(31) 



4-> <— ^ 

where Z, = Zt -h Z,, and 

^ In (31), the reflection coefficients are conveniently those at the endpoints (note: Z{zr) =ZtR'> 
Z{zl) =ZrL), and the entire expression may be written out explicitly in terms of the travehng wave 
solutions of the homogeneous equation (23). It is noted in this connection that one may put 
IZ = I]{, II = 1^, since only normalized quantities appear in the formula. An analogous result for 
the Green's function Z(z, z') follows from duality considerations. 

Equation (31) simplifies considerably for certain special cases. In an infinitely extended 

region, r(z/J =r(z/;) =0, and one finds 

r (., .')- ^^' ^^<^ '^~ ^^' ^ ^ ^^' ^^"^ '^' ^'' ^ ^ (32) 

K(z') 
which formula is also deduced directly from (26). In a semi-infinite region (with zl^ — ^, i.e. 

nz,)=o), 

Ziiz')I-{z')hi'') 
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Returning now to the problem at hand after these general considerations, we see from (10) or 
(11) that the homogeneous equations for Y> 1 are 



|^-hA:2 [1- J'^'_ ]}r^0 for the E modes, 



(34a) 



^ + A:2 [1-M] j r = for the H modes, 



(34b) 



which results follow also from (23) and its dual, in view of the definitions of k, Zc in (6). The travel- 
ing wave solutions for the H modes evidently comprise the exponentials exp [± i VA:^ — kt^z] , 
independent of the medium constants, whereas their form for the E modes depend on the details 
of the stratification X(z). 



3. Application to a Semi-Infinite Medium With Linearly Varying Electron 

Density 

A simple illustration of the preceding results is now given for the case of a semi-infinite, nearly 
uniaxial medium with a linear variation of X with z as shown in figure 2. In a subsequent paper, 
it is planned to cover further details of the derivation, together with applications to various other 
nonuniform configurations. 

The source is taken as the transverse electric current element located at 2 = 2:' in the vacuum 
half space z < zt, = — — (see sec. 2.3). The independent functions /+ and /~ for both E and H modes 
are simple exponentials in the vacuum region (the pertinent equation is the same as (34b) since 
X = 0). In the linear half space the H mode functions remain unaltered; however, the E mode 
differential equation is solved in terms of Whittaker functions. Throughout the entire space, 
the H mode Green's function is 



Z"{z,z')- 



CO^to 



2Vk^-kf 



^iV¥^\z-z' I 



(35) 



a resuk obtained at once from the dual of (32), with Yi'iz') =2F^' = 2/c7aj/xo [see (6)]. 

The E mode Green's function, constructed from (33) (with zr = Zb), may be shown to be given by 



weo 



Y{z,z')-- 



2VF-k'j 
for z < Zh, and by 



|-^,-vPr7^^_^'l_^ 



iVk^'-kfW,{zb)-JVl{zb) 
iV¥^Wx{zb)^W((zb] 



. ^ - i\k2 - mz + 2' - 2z ) "I 



Y{z, z')= icoeoe-^^^'-^y -h^ 



WAz) 



TRANSVERSE 
BOUNDARY PLANE 



POINT SOURCE ^ 

^^ 



iVk^-kfWi{zb)'hW[{zb) 



I I SLOPE A 

-^ I X = 1+^2 

I t 

FREE SPACE H X | 




(36a) 



(36b) 



Figure 2. Semi-infinite linear model. 
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for z > Zii, where Zi, locates the boundary plane separating the two semi-infinite regions shown 
in figure 2. JVi and Wz are hnearly independent Whittaker functions defined by [Slater, 1960] 



W,{z) = WjI_ 1 {-i2kz) 
2Ji ' ' 

W.i{z) = r-«| 1 {i2kz) 
2j.V 



> 



where a- is the derivative of X with respect to z and the primes on W\ in (36a, b) denote 



dzi, ' 



(37) 



The 



reflected wave in (36a) appears to originate from the image point at 2zh — z', and tlie factor muhi- 
plying the second exponential is the reflection coefficient T{zb), which was derived from a knowl- 
edge of Z{zb)\ i.e., of the voltage and current solutions in the plasma half-space. Since W\{z) 
may be shown to satisfy the radiation condition at z ^ oo, and with Z, {zu) — V A:^ — kf I ojeo in vacuum, , 
the form of F (zb) in (36a) follows. The formula in (36b) is obtained at once from the continuity 
of Y(z, z') at zi) and the above-mentioned behavior of Wiiz) at z^ oo. 

The uniaxial portion of the far field solution in the vacuum half-space, obtained from an asymp- 
totic evaluation, is found to be given by [Gross, 1964] . 



Eou 



■—. — cos d cosi//(l — 1 ) 

477 r 



ik . , e'^"- 
477 r 

H,\,u —~1 — cos cos i//(l — F')- — 



Eru — H,u — 



(38) 



where (r, 8, i//) are the spherical coordinates of the observation point relative to the source; Eou, 
E^w, etc., are the spherical field components, and F' is the reflection coefficient 



^, _ ik cos W,{Zb) - W[{Zb) ^_,-..,.._.^^. eos 

ik COS Wx{zb) -h W[{zu) 



(39) 



The angle B here is measured from the negative z-axis. 

The first order gyrotropic corrections to the far field are found from section 2.3b and turn out 
to be [Gross, 1964]: 



EoG — 

17 

^il/G — ' 

HdG = 



— sm lb / 

477 ^ r 



oikr 



. COS 6 cos \1) .y 

477 ^ r 

— iJ^ _ ^ikr 

- — cos 6 cos \b '/ 

477 ^ r 



ik 



oikr 



H4,G = - — sin i// 7 

477 r 



ErG — HrG — 



(40) 
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where 



'f=- 



^-i2k{z'-Zf)) cos d 



YW,{zb) + ik cos 9 Wiizb). 



("dzV- 

JZh 



.z")e' 



/k'-ky-zi) 



w[iz'r 



(41) 



To account for the modification due to (40), one may define a new refliection coefficient FJ, 



r;=r'-ysec ^tam// 



(42) 



which highhghts the skewing of the far field pattern with i//. The power patterns are similarly 
modified by the gyrotropic correction [Gross, 1964]. 

The integral in (41) may be evaluated approximately with the assumption of slight losses, 
thereby effectively truncating the integration interval. When the rate of variation of X with z is 
large, it is found that 



'/=- 



A:F(l-hcos 6)^ 



^-i2k{z' -Zfy) cos 6 



(43) 



from which 



EoG — 



ico/xo e" 



sin i// 



47r r (l-hcos^)3 



^-i2kiz'-Zfy) cos t 



(44) 



and similarly for other components. Details of this calculation, as well as of the more involved 
small slope approximation, will be presented in a subsequent paper. Although a large slope is 
admitted, the inequality ^ < < kY must be satisfied in order to retain the perturbation character 
of the gyrotropic correction. Graphs of the angle dependent parameters of Eeu, Eec and E^g are 
shown in figures 3 and 4. 
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Figure 3. Uniaxial field component for the semi-infinite 
linear model. 
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The ordinate represents {Eeul in the xz plane muhiplied by the factor , [ | |_j j^g abscissa 0. 
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Figure 4. Relative gyrotropic corrections to field 
components. 

The ordinate represents the pertinent field quantity, lEecl or lE^id, multipHed by 
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4. Discussion 

Figure 3 contains graphs descriptive of \E0u\ versus in the xz plane (i// = 0°) for various values 
of the normalized slope parameter ^X, where X is the free-space wavelength. The curves repre- 
sent the angular function of Esu in (38), with the source taken at the boundary {j.z' =— 1). Values 
of ^ X were selected to demonstrate different characteristics. The ^A ~ curve corresponds 
to the envelopes of the peaks of the multilobe pattern obtained for small values of ^k. The 
^X = 7r/5 curve exhibits a four lobe pattern. This curve was plotted by using an approximation 
for F' valid when the rate of density variation in the medium is small, and contains nulls that go 
to zero. For comparison, a better approximation for the region near 6^ = is shown by a dashed 
curve. The ^k = 7r curve was plotted using exact values of the Whittaker functions in the ex- 
pression for F'. The curve has two lobes and a broader null that does not dip to zero. As ^X 
increases further the lobes disappear and the pattern approaches that shown for ^X^oo. The 
free space curve is also shown for comparison with the ^X ^ oo curve. 

The multilobe patterns which are found for small values of the slope parameter are typical 
of highly reflective media. The pattern for an infinite value of the parameter is similar to the free- 
space curve shown in the figure, and is typical of a weakly reflecting region. Even when the 
source is moved away from the boundary these characteristics are basically unaltered [Gross, 
1964]. The curves show, therefore, that more reflection is obtained for the small slope case, 
whereas regions with large slope reflect only a small amount of energy; the rest penetrates deeper 
into the medium. These results have also been examined in detail from a physical standpoint, 
for example, using ray theory [Gross, 1964], and have been found consistent with such viewpoints. 

Figure 4 illustrates the variation of the angular factors of lEea] and \E^},a\ with 6 in their plane 
of maximum magnitude, i// = 90°. These curves are plotted for the large slope approximation and 
are independent of the source point z\ since z' enters only in an imaginary exponential factor. The 
phase varies with z\ producing interference or enhancement of the uniaxial field components. 
\E0a\ reaches a maximum at 6^ = 90°, whereas, l^^r/l is almost independent of the angle for 6 < 80° 
and decreases to zero beyond. The magnitudes shown in the figure are only relative, since the 
absolute numerical contribution depends on the value of F, which is considered suitably large for 
the purpose here. 

Although lEecl reaches a maximum at 6^ = 90°, it can be shown [Gross, 1964] that the first- 
order gyrotropic power is zero at this angle. In fact the Poynting vector is zero throughout the 
xz and yz planes where the only contributions from the gyrotropic corrections must be from higher 
order terms. It can also be shown that the gyrotropic power peaks in the i// = 45°, 135° planes 
and is independent of z' . With respect to 6, the power is found to peak at 6 = IT for large values 
of ^ X. These results indicate the skewness of the gyrotropic corrections. 
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